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Abstract 


Genomic imprinting is found in marsupial and eutherian mammals, but not in monotremes. While the primary regulator of 
genomic imprinting in eutherians is differential DNA methylation between parental alleles, conserved imprinted genes in 
marsupials tend to lack DNA methylation at their promoters. DNA methylation at eutherian imprinted genes is mainly cata- 
lyzed by a DNA methyltransferase (DNMT) enzyme, DNMT3A. There are two isoforms of eutherian DNMT3A: DNMT3A and 
DNMT3A2. DNMT3A2 is the primary isoform for establishing DNA methylation at eutherian imprinted genes and is essential 
for eutherian genomic imprinting. In this study, we investigated whether DNMT3A2 is also present in the two other mam- 
malian lineages, marsupials and monotremes. We identified DNMT3A2 in both marsupials and monotremes, although im- 
printing has not been identified in monotremes. By analyzing genomic sequences and transcriptome data across vertebrates, 
we concluded that the evolution of DVMT3A2 occurred in the common ancestor of mammals. In addition, DNMT3A/3A2 
gene and protein expression during gametogenesis showed distinct sexual dimorphisms in a marsupial, the tammar wallaby, 
and this pattern coincided with the sex-specific DNA methylation reprogramming in this species as it does in mice. Our results 
show that DNMT3A2 is present in all mammalian groups and suggests that the basic DNMT3A/3A2-based DNA methylation 
mechanism is conserved at least in therian mammals. 
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Significance 


While eutherian genomic imprinting is mainly regulated by DNA methylation, marsupial imprinted genes tend to lack 
promoter DNA methylation. In this study, we asked whether this difference is caused by the absence of a primary en- 
zyme, DNMT3A2. We confirmed that DNMT3A2 is present in all mammalian lineages. As the primary enzyme must have 
been evolved in the common ancestor of mammals, our result suggests that the difference in DNA methylation levels 
between marsupial and eutherian imprinting is caused by the associated proteins rather than the enzyme itself. 


Introduction 


Genomic imprinting is an epigenetic process that causes 
differential expression of a subset of genes depending on 
their parental origin. Among vertebrates, genomic imprint- 
ing has been found in therian mammals, namely eutherians 
and marsupials, but there is no evidence for imprinting in 


monotremes or nonmammalian vertebrates (Pask et al. 
2009; Renfree et al. 2009; Frésard et al. 2014; Griffith 
et al. 2016). Mammalian genomic imprinting is therefore 
thought to have evolved after the therian-monotreme di- 
vergence ~184 million years ago (Cuneo et al. 2013). 
Genomic imprinting is known to regulate mammalian de- 
velopment and maternal behavior (Lefebvre et al. 1998; 
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Reik and Walter 2001; Plagge et al. 2004; Keverne and 
Curley 2008; Renfree et al. 2009; Ferguson-Smith 2011), 
and disruption of imprinting causes many developmental 
defects such as placental abnormalities (Ono et al. 2006; 
Frost and Moore 2010; Piedrahita 2011). Although a num- 
ber of hypotheses have been proposed for the roles of im- 
printed genes, questions remain as to how and why 
genomic imprinting evolved in therian mammals. As this 
phenomenon must have evolved in the common ancestor 
of therian mammals, characterizing similarities and differ- 
ences in imprinting mechanisms between marsupials and 
eutherians, and comparing these to monotremes allow us 
to begin to characterize the ancestral mechanisms and to 
define how mammalian imprinting evolved. 

The majority of the known-eutherian imprinted genes 
have certain regions marked by DNA methylation that differs 
in its methylation status between maternal and paternal al- 
leles. Such regions are known as differentially methylated re- 
gions (DMRs; Bartolomei and Ferguson-Smith 2011). 
Differential DNA methylation status between the parental al- 
leles is considered as the primary mechanism for genomic 
imprinting in eutherians. However, only three DMRs have 
so far been reported in marsupials (Weidman et al. 2006; 
Suzuki et al. 2007, 2018; Smits et al. 2008). Four conserved 
imprinted genes, MEST, IGF2R, L3MBTL, and HTR2A, which 
have a promoter DMR in eutherians (Stéger et al. 1993; 
Lefebvre et al. 1997; Riesewijk et al. 1997; Li et al. 2002, 
2004; Joshi et al. 2016), lack a DMR on their promoter region 
in marsupials (Suzuki et al. 2005, 2018; Das et al. 2012). It is 
currently unknown how their imprinting status is established 
and maintained without a promoter DMR. In addition, mar- 
supial imprinting is generally considered less stable than eu- 
therian imprinting as it shows leaky expression from the 
silenced allele (Renfree et al. 2008). This indicates that eu- 
therian mammals acquired additional molecular machineries 
for establishing more stable inactive marks like DNA methy- 
lation on promoters of genes such as MEST and /GF2R after 
the marsupial-eutherian split 166 million years ago. 

In eutherians, DNA methylation is catalyzed by the DNA 
methyltransferase (DNMT) family of enzymes including 
DNMT3A, DNMT3B, and DNMT3L (Auclair and Weber 
2012; Saitou et al. 2012; Sadakierska-Chudy et al. 2015). 
To ensure a sex-specific DNA methylation pattern is trans- 
mitted in gametes, the genome-wide DNA methylation sta- 
tus is reprogrammed in developing germ cells in both 
marsupial and eutherian mammals (Seisenberger et al. 
2012; Monk 2015; Ishihara et al. 2019). In mice, this repro- 
gramming is achieved by loss of parental DNA methylation 
marks followed by the resetting of DNA methylation by 
Dnmt3a, Dnmt3b and Dnmt3l in a sex-specific manner 
(Okano et al. 1999; Kaneda et al. 2004). In female mice, 
conditional KO of Dnmt3a reduces DNA methylation at 
the DMR of the conserved paternally expressed gene 
Mest (Kaneda et al. 2004), indicating that the DMR of 


this gene is established by the action of DNMT3A in euther- 
ians. As MEST in marsupials is imprinted but does not have 
a DMR, comparisons of DNMT3A-based mechanisms 
between marsupials and eutherians could elucidate the dif- 
ferences in DMR-based imprinting between the two mam- 
malian lineages. 

There is a novel short isoform of DNMT3A, DNMT3A2, 
identified in eutherians (Chen et al. 2002; Sakai et al. 
2004). Further characterization of DNMT3A2 demon- 
strated that it is a more effective isoform of DNMT3A at es- 
tablishing DNA methylation on imprinted genes (Sakai et al. 
2004; Nimura et al. 2006; Suetake et al. 2006; Ma et al. 
2015; Thakur et al. 2016). Marsupials also have DNMT3A, 
DNMT3B and DNMT3BL (Yokomine et al. 2006; Ishihara 
et al. 2019), but the presence of DNMT3A2 is unknown. 
It is therefore possible that the evolution of DNMT3A2 oc- 
curred in the ancestor of eutherian mammals and caused 
the differences observed in the DNA methylation status of 
imprinted genes between eutherians and marsupials. 

In this study, the evolution DNMT3A2 was investigated 
through the identification of isoforms of DNMT3A across 
vertebrates focusing on the two other mammalian lineages, 
marsupials and monotremes. By analyzing transcriptome 
data of a marsupial (the koala: Phascolarctos cinereus), a 
monotreme (platypus: Ornithorhynchus anatinus), and 
non-mammals (chicken: Gallus gallus; the bearded dragon: 
Pogona vitticeps), we evaluated the presence of DNMT3A2 
across vertebrates. We further compared the gene and pro- 
tein expression patterns of DNMT3As in the developing go- 
nads of the marsupial tammar wallaby (Macropus eugenii) 
with the existing DNA methylation data for this species to 
determine conserved function. Here we report that the 
DNMT3A2 occurred in all three mammalian groups so it 
must have evolved in the common ancestor of mammals 
and that the basic DNMT3A/3A2-based mechanism ap- 
pears to be conserved at least in therian mammals. 


Results 


Tammar Wallaby has the Short Form of DNMT3A, 
DNMT3A2 


The 2,865 bp of tammar DNMT3A candidate sequence 
was identified by a blast search of the tammar genome 
database (Wallabase: https:/wallabase.org). The com- 
parison of the putative tammar DNMT3A with mouse 
Dnmt3a/Dnmt3a2 identified highly conserved protein 
coding exons (fig. 1A and supplementary fig. $1). 5’-rapid 
amplification of CDNA end (RACE) reactions using a specific 
primer against the highly homologous region (fig. 1A) were 
then performed with adult tammar testis to characterize 
the marsupial DNMT3A isoforms. The 5’-RACE products 
showed two distinct bands (fig. 18), indicating that there 
were two isoforms. After cloning and sequencing of each 
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Fic. 1.—ldentification of the tammar DNMT3A2. RACE experiments identified two isoforms of DNMT3A including DNMT3A2 in the tammar. (A) Exon 
structures of mouse Dnmt3a/3a2 and putative tammar DNMT3A. Conserved exons used for primer design were represented by green-colored boxes. 
(B) Identification of TSS of the tammar DNMTS3A isoforms. 5’-RACE and nested RACE primers were represented by the black arrows. Asterisks represent 
RACE products coding partial tammar DNIVT3A sequences. Boxes represent exons identified by sequencing of the RACE products. Orange-colored box shows 
DNMT3A2-specific exon. (C) Identification of the full length of the tammar DNMT3A and DNMT3A2. Standard PCR primers designed for characterizing each 


isoform and 3’-RACE primers were represented by the black arrows. 


RACE product, two transcript variants of DNMT3A were iden- 
tified. The first exon of the shorter isoform was located be- 
tween exon 6 and exon 7 of the longer isoform. Based on 
the exon structure, we renamed the shorter isoform as 
DNMT3A2 (fig. 1C). Combining specific primers against the 
DNMT3A2 and 3’-RACE reactions confirmed that the tammar 
DNMT3A2 is also a protein coding transcript (fig. 1). 


Platypus also has DNMT3A2 and a Unique Isoform of 
DNMT3A, DNMT3A0a 


The 2,865 bp of platypus DNMT3A candidate sequence was 
identified by a blast search of the new platypus genome 
(mOrnAna1.p.v4). 5’-RACE reactions were then performed 
using adult platypus testis to identify monotreme DNMT3A 
isoforms (fig. 2A). The 5’-RACE products showed four ampli- 
cons, indicating that there were four isoforms. After sequen- 
cing each RACE product, four transcript variants of DVMT3A 
were identified. Isoform 4 encoded the shortest form of 
DNMT3A and its first exon was located between exon 6 
and exon 7 of the longest form. Based on the exon structure, 


we renamed the shortest isoform as DNMT3A2 (fig. 28). 
Using specific primers against the DNMT3A2 and isoform3, 
3’-RACE reactions confirmed that the platypus DNMT3A2 
and isoform3 were protein coding transcripts (fig. 2C). 
Isoform 3 contained a unique exon and encoded a unique 
protein. A BLAST search using NCBI BLAST did not reveal 
any similar sequences with this unique exon. Thus, we 
named isoform 3 as a novel isoform, DNMT3A O. anatinus 
specific isoform, DNMT3AOQa (fig. 2C). 


The Transcription Start Site of DNMT3A2 is Conserved 
Across Mammals 


To identify conserved cis-regulatory elements at the tran- 
scription start site (TSS) of DNMT3A2 in mammals, genomic 
alignment was performed and visualized with AliTV pro- 
gram (Ankenbrand et al. 2017). The TSS of DNMT3A2 
was highly conserved across mammals (fig. 3). Note that 
the TSS of platypus DNMT3A2 also has high homology 
with a region in the chicken DNMT3A locus. However, 
the range of the region is <62 bp and the sequence does 
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represent race products coding partial platypus DNMT3A sequences. Boxes represent exons identified by sequencing of the RACE products. Orange- and 
red-colored boxes show DNMT3A2-specific and DNMT3AOa-specific exons, respectively (C) Identification of the full length of the platypus DNMT3As. 
Standard PCR primers designed for characterizing each isoform and 3’-RACE primers were represented by the black arrows. 


not encode the first exon of mammalian DNMT3A2. 
The TSS of mammalian DNMT3A2 showed no homology 
with the DNMT3A locus of the bearded dragon lizard 
(P. vitticeps). 

To identify common cis motifs around the TSS of 
DNMT3A2, the highly homologous region identified across 
mammals was subjected to multiple alignments by 
clustalW. As a result, highly conserved sequences around 
the TSS of DNMT3A2 across mammals were identified 
(fig. 4A; asterisks). Regions containing five or more highly 
homologous sequences were compared with known tran- 
scription factor binding motifs using the TOMTOM pro- 
gram in the MEME suite. From this analysis, two common 
regions were found that had significant matches with 


F-values <0.05 and q-values <0.05 (fig. 4A; dotted boxes). 
Region 1 contained upstream transcription factor 2 (USF2) 
and USF1 binding motifs (fig. 48), and region 2 contained a 
T-box transcription factors (TBX15, TBX1, TBX4, TBX5) 
binding site and a Max dimerization protein (MGA) binding 
site (fig. 4C). Region 2 was downstream of the TSS and 
within the first exon of mammalian DNMT3A2. 


DNMT3A2 is a Mammal Specific Isoform of DNMT3A 


Published transcriptome data sets of gonads derived from 
chicken and the bearded dragon were analyzed to identify 
potential isoforms by visualizing splicing junctions. Koala 
and platypus transcriptome data sets were also analyzed 
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and DNMT3A2 are shown as black box connected with black lines in the 
aqua boxes. 


to evaluate the isoform detection method based on visual- 
izing splicing junctions. Visualizing splicing junctions suc- 
cessfully identified the short form DNMT3A2 from both 
the koala and platypus transcriptome data (fig. 5A and 8). 
In the platypus, the presence of DNMT3AOa was also con- 
firmed in both male and female gonadal transcriptome 
data by detecting reads mapped to the unique exon of 
the isoform (fig. 5B). However, chicken and bearded dra- 
gon data sets clearly lacked the short form (fig. 5C and D). 


Expression Patterns of DNMT3A/A2 During 
Gametogenesis are Conserved Across Therian Mammals 


Real-time reverse transcription-polymerase chain reaction 
(RT-PCR) experiments using developing gonads from the 
tammar were performed to confirm the correlation be- 
tween DNMT3A/A2 expression and the existing data on 
tammar DNA methylation reprogramming during gameto- 
genesis. During oogenesis, tammar DNMT3A expression 
was upregulated at day 200 pp, after folliculogenesis be- 
gins (Mattiske et al. 2002; fig. 6A). During spermatogen- 
esis, tammar DNMT3A_ expression was_ significantly 
upregulated during mitotic arrest when most of the DNA 
methylation is re-established (Ishihara et al. 2019) (fig. 
6B). In both sexes, tammar DNMT3A2 showed significantly 
higher expression at day 10 pp compared with the other 
ages examined in this study (P<0.05). After day 25 pp, 
tammar DNMT3A2 expression was downregulated in 
both sexes (fig. 6C and D). In females, tammar DNMT3A2 


expression was significantly upregulated (P<0.05) at day 
120 pp when DNA remethylation begins (Ishihara et al. 
2019) and then gradually downregulated as development 
proceeds (fig. 6C). In contrast to females, male tammar 
DNMT3A2 expression was constant during prespermato- 
genesis and significantly decreased (P< 0.05) in the adult 
testis (fig. 6D). 


DNMT3A/3A2 Protein became Nuclear at the Time of 
DNA Remethylation in Both Sexes 


To further confirm the correlation between DNA methyla- 
tion reprogramming and DNMT3A/3A2 expression, immu- 
nohistochemistry (IHC) was performed. Based on our 
western blotting analysis, the antibody used in this study 
could detect both DNMT3A = and = DNMT3A2 
(supplementary fig. S2). In developing gonads, tammar 
DNMT3A/3A2_ protein localization was always nuclear 
(fig. 7). Within ovaries, at day 16 pp, DNMT3A/3A2 was re- 
stricted to somatic cells (fig. 7A). By day 60 pp, when early 
folliculogenesis has commenced, several oocytes within 
primordial follicles adjacent to the cortex were DNMT3A/ 
3A2-positive (fig. 7B). At day 120 pp, when DNA remethy- 
lation begins (Ishihara et al. 2019), some oocytes contained 
nuclear DNMT3A/3A2 but others were negative, seemingly 
with no respect to oocyte size or location (fig. 7C). At day 
230 pp, some smaller oocytes within primary follicles had 
nuclear DNMT3A, as did most larger oocytes within sec- 
ondary and tertiary follicles (fig. 7D). 

At day 26 pp, the majority of germ cells were DNMT3A/ 
3A2 negative (fig. 7£). At day 60 pp when DNA remethy- 
lation has commenced (Ishihara et al. 2019), it appeared 
there were more DNMT3A/3A2-positive cells compared 
with day 26 pp (fig. 7F). Day 110 pp testes contained simi- 
lar numbers of DNMT3A/3A2-positive cells to day 60 pp 
(fig. 7G). However, in day 300 pp testes, DNMT3A/ 
3A2-positive germ cells were not evident compared with 
day 110 pp (fig. 7H). 


Discussion 


The short form of DNMT3A, DNMT3A2 is critical in estab- 
lishing DMR-based imprinting in eutherians. As most mar- 
supial imprinted genes so far identified do not have 
DMRs, we predicted that the evolution of DNMT3A2 oc- 
curred after the eutherian-marsupial split 166 million years 
ago. However, DNMT3A2 was present not only in marsu- 
pials but also in monotremes, in which genomic imprinting 
has not been identified. No DNMT3A2 was identified in the 
chicken or the bearded dragon, but the TSS of DNWMT3A2 
was conserved across all mammals, suggesting that the ac- 
quisition of the new promoter for DNMT3A2 expression oc- 
curred in the common ancestor of mammals before the 
monotreme-therian split over 184 million years ago 
(Cuneo et al. 2013; fig. 8). 
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common cis-acting elements. (A) Multiple alignments of the highly conserved region around the TSS of DNMT3A2 in mammals. The first exon of DNMT3A2 is 
shown in red-colored letters. There are two regions (regions 1 and 2) which had significant matches (E-value <0.05 and q-value <0.05) based on TOMTOM 
program in MEME suite. (8) Motif structure of the region 1. The region contained upstream transcription factor 2 (USF2) and USF1 binding motifs. (C) Motif 
structure of the region 2. The region contained T-box transcription factors (TBX15, TBX1, TBX4, TBX5) binding site and Max dimerization protein (MGA) bind- 


ing site. 
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Fic. 5.—Visualizations of splicing junctions of DNMT3A confirmed DNMT3A2 in koala and platypus but not in other vertebrates. Visualizations of splicing 
junctions of vertebrate DNMT3A identified DNMT3A2 in koala and platypus but not chicken and bearded dragon. Published gonadal transcriptome data sets 
derived from koala, platypus, chicken, and a dragon lizard were mapped against their genomic sequences and splicing junctions of DNMT3A were visualized 
by sashimi plot. Splicing junctions were shown as lines connecting two exons and numbers of junction reads were shown on the junction lines. Exon structures 
of DNMT3A were shown as black boxes under the sashimi plot. (A) Koala. In testis, koala transcriptome data analysis successfully detected the short form, 
DNMT3A2 (orange circle). (B) Platypus. In both sexes, platypus transcriptome data analysis successfully detected the short form, DNMT3A2 (orange circle). 
(C) Chicken. Chicken transcriptome data analysis could not identify DNMT3A2 in gonads. (D) The central-bearded dragon, Pogona vitticeps. 


Transcriptome data analysis did not identify DNMT3A2 in gonads. 
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Fic. 6.—Expression of DNMT3A and DNMT3A2 during tammar gametogenesis. DNMIT3A/3A2 gene expression during gametogenesis was consistent 
with DNA methylation reprogramming in the tammar. (A) DNMT3A expression during spermatogenesis. (B) DNMT3A expression during oogenesis. 
(C) DNMT3A2 expression during spermatogenesis. (D) DNMT3A2 expression during oogenesis. Relative gene expression against two housekeeping genes, 
HMBS and TBP, was shown as bar graph. Error bar represents standard error of the mean (SEM). Bars labelled with the same letters represent values that do not 
differ significantly by Tukey-Kramer's multiple comparison test at P< 0.05. Sample numbers of each stage are represented in parentheses. The corresponding 
developmental stages of male and female tammar germ cells based on pervious publications (Ullmann et al. 1997; Ishihara et al. 2021) are shown at the 
bottom. DNA methylation data were based on a previous study (Ishihara et al. 2019). 


The highly homologous region around the TSS of 
DNMT3A2 across mammals contained a number of consen- 
sus binding motifs. The putative USF1/2 binding sites were 
located in the TSS of the tammar DNMT3A2 and upstream 
of the TSS of DNMT3A2 in the other mammals examined in 
this study. USFs generally bind to genomic regions close to 
TSS of target genes and regulate their expression 
(Rada-lglesias et al. 2008). For example, USF1 binds to the 
first exon of the N4bp2/71 gene and regulate its expression 
in mice (Watanabe et al. 2019). This suggests that USF1/2 
binding may contribute to the expression of DNMT3A2 in 


mammals. Moreover, we found a putative binding site of 
T-box transcription factors such as TBX15 and TBX5. The 
putative TBX factor binding site was downstream of the 
TSS of DNMT3A2. However, Tbx5 binds to a motif located 
downstream of the TSS of Fgf70 and controls its expression 
in developing mouse limbs (Delgado et al. 2021). Thus, 
T-box transcription factors regulate target gene expression 
even when they bind to binding sites downstream of the 
TSS. These data suggest that these T-box transcription fac- 
tors may regulate DNMT3A2 expression in mammals as 
seen in mouse Fgf70. TBX15 is known to be required for 
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Fic. 7—DNMT3AV3A2 protein localization during tammar gametogenesis. Immunohistochemistry of DNMT3A/3A2 during oogenesis and spermatogen- 
esis revealed sex-specific DNMT3AV/3A2 upregulation in developing germ cells. (A) Ovary at day 16 pp. (B) Ovary at day 60 pp. (C) Ovary at day 110 pp. 
(D) Ovary at day 192 pp. (E) Testis at day 26 pp. (F) Testis at day 60 pp. (G) Testis at day 120 pp. (H) Testis at day 300 pp. Black and red arrows represent 
germ cell with DNMT3A/3A2 staining and DNMT3AV/3A2 negative germ cells, respectively. Brown color: positive staining. Scale bars: 20 um. 


skeletal development of the limb, vertebral column, and 
head (Singh et al. 2005). In mice, Dnmt3a/3a2 deficiency 
shows defects in osteoclast differentiation (Nishikawa 
et al. 2015). The other T-box transcription factor, TBX5, is 
involved in cardiovascular development (Steimle and 
Moskowitz 2017), and inactivating DNMT3A/3A2 contri- 
butes to cardiovascular disease (Sano et al. 2018). 
Therefore, the short isoform, DNMT3A2, in addition to es- 
tablishing DNA methylation-based imprinting in germ cells, 
may have a functional role in mammalian bone and cardio- 
vascular development. Further characterization of the 


molecular mechanisms that regulate DNMT3A2 expression 
could elucidate possible links between DNMT3A2, T-box 
transcription factors and mammalian development. 

A cofactor of DNMT3A, DNMT3L is essential for adding 
DNA methylation on selected gene loci, including imprinted 
genes in eutherians (Bourc’his et al. 2001; Hata et al. 2002; 
la Salle et al. 2007; Glass et al. 2009). Across mammals, 
DNMTBL is conserved in both marsupials and eutherians 
but not in monotremes (Yokomine et al. 2006; Renfree 
et al. 2009). These data suggest that monotremes lack an 
ability to recruit DNA methylation on selected gene loci to 
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Fic. 8.—Evolution of DNMT3A2 occurred in the common ancestor of 
mammals. Phylogenetic tree and DNMT3s. DNMT3 family member pro- 
teins of each animal are listed under the silhouette of each animal. The im- 
portant isoform of DNMT3A, DNMT3A2 has evolved in the common 
ancestor of mammals (red letters). The platypus lacks a cofactor of 
DNMT3A, DNMT3L, but the platypus has the platypus-specific isoform of 
DNMT3A, DNMT3A0Oa (blue letters). 


establish genomic imprinting. However, in this study, the 
novel platypus-specific isoform DNMT3A0Oa was identified. 
Although the exact function of DNMT3A0Oa is currently un- 
known, this novel isoform may increase the capacity for the 
platypus to add DNA methylation onto their genome. The 
capacity of DNMT3AOa to add DNA methylation on platy- 
pus genome needs to be determined. 

In mice, Dnmt3a/3a2 expression shows distinct sexual di- 
morphisms during gametogenesis and correlates with DNA 
methylation reprogramming in both sexes (la Salle et al. 
2004; Lucifero et al. 2007). Our analysis showed that pat- 
terns of tammar DNMT3A and DNMT3A2 expression 
were also different between the sexes during gametogen- 
esis and correlated with the known DNA methylation pat- 
terns across gametogenesis in the tammar. Although 
the previous mouse studies did not separate Dnmt3a2 
expression from Dnmt3a expression in developing germ 
cells (la Salle et al. 2004; Lucifero et al. 2007), the combined 
expression of mouse Dnmt3a/Dnmt3a2 and tammar 
DNMT3A/DNMT3A2 have similar patterns during gameto- 
genesis. Our IHC analysis further confirmed that DNMT3A/ 
3A2 localization also showed sex-specific differences 
consistent with DNA methylation reprogramming in this 
species. These data suggest that the sex-specific global 
DNA methylation reprogramming in therian mammals 
is mediated by a conserved DNMT3A/DNMT3A2-based 
mechanism. Although we could not confirm DNMT3A/ 
3A2 expression in developing platypus gonads due to lim- 
ited access to samples, it would be interesting to determine 
if this pattern is conserved across all mammals. 

As well as DNMT3A/3A2 expression, the expression pat- 
tern of its co-factor DNMT3L in the tammar is also similar to 


mouse Dnmts3! (Ishihara et al. 2019). Therefore, based on 
this study and previous studies (Yokomine et al. 2006; 
Renfree et al. 2009; Ishihara et al. 2019), we confirmed 
that marsupials have all the known molecular machineries, 
DNMT3A/3A2 and DNMT3BL, required for catalyzing DNA 
methylation on imprinted genes. Despite this, some marsu- 
pial imprinted genes lack promoter DNA methylation. For 
example, in the tammar, the conserved imprinted gene 
MEST lacks a DMR on its promoter (Suzuki et al. 2005, 
2018; Das et al. 2012). However, in mice, Mest DMR estab- 
lishment in oocytes occurs in the later stages of oogenesis 
(Lucifero et al. 2004) and during tammar oogenesis, late 
upregulation of DNMT3A relative to DNMT3A2 was ob- 
served. As the combined expression of mouse Dnmt3a/ 
Dnmt3a2 and tammar DNMT3A/DNMT3AZ2 has similar pat- 
terns during oogenesis, this suggests that the lack of DNA 
methylation on the promoter of marsupial MEST might be 
caused by differences in DNMT3A and/or its associated pro- 
teins between marsupials and eutherians. However, con- 
firmation of this requires examination of the specific 
expression pattern of Dnmt3a relative to Dnmt3a2 in the 
mouse during oogenesis. 

This study confirms that the overall mechanism and ex- 
pression pattern of DNMT3A/DNMT3AZ2 is conserved 
among therian mammals. However, given this conserva- 
tion, it remains unclear as to why some marsupial imprinted 
genes lack DMRs on their promoters. DNMT3A2 was also 
present in monotremes. In addition, a novel DNMT3A iso- 
form was identified in monotremes which may contribute 
to establish the global DNA methylation levels. Further 
examination of monotremes could provide important in- 
sights into evolution of mammalian genomic imprinting. 


Materials and Methods 


Animals 


Tammar wallaby (M. eugenii) samples of Kangaroo Island, 
South Australia origin, were collected from captive animals 
from our colony maintained by the University of 
Melbourne. Gonads were snap-frozen immediately after 
dissection. Platypus (O. anatinus) tissues were collected 
from wild-caught animals under permits from NSW Parks 
and Wildlife and ethically approved by the University of 
Melbourne Animal Ethics committees. All animal handling 
and husbandry were in accordance with the National 
Health and Medical Research Council of Australia (2013) 
guidelines. 


Gonadal RNA Extraction and cDNA Synthesis 


Snap-frozen gonads were used for RNA extraction using 
the GenElute Mammalian total RNA Miniprep. Kit 
(Sigma-Aldrich, St Louis, MO, USA) following the manu- 
facturer’s instructions. The extracted RNA was treated 
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with the DNA-free DNase treatment and removal kit 
(Thermo Fisher Scientific, Waltham, MA, USA) to remove 
residual genomic DNA then 200 ng of RNA was used as a 
template for cDNA synthesis using SuperScript IV First 
strand Synthesis System (Invitrogen, Carlsbad, CA, 
USA). Approximately 1 ug of RNA was used as a template 
for cDNA synthesis using a SMARTer RACE 5’/3’-kit 
(Clontech, Mountain View, CA, USA). 


5’- and 3’-Rapid Amplification of CDNA Ends 


To determine the TSSs of the marsupial and monotreme 
DNMT3A transcripts, 5’-RACE experiments were per- 
formed using the SMARTer RACE 5//3’-kit (Clontech). 
Putative tammar and platypus DNMT3A sequence were 
searched by a blast search of the wallaby genome 
(Wallabase: https:/wallabase.org) and the platypus gen- 
ome (mOrnAna1.p.v4), respectively. After identifying the 
putative sequences, gene-specific primers were designed 
for the downstream RACE experiments. The first round 
RACE reaction was performed with adult tammar testis or 
adult platypus testis cDNA using SeqAmp DNA 
Polymerase (Clontech) with the gene-specific primers 
(supplementary table S1). The nested 5’- and 3’-RACE 
was performed by GoTaq DNA polymerase (Promega, 
Fitchburg, WI, USA) and the RACE products cloned using 
pGEM-T Easy Vector (Promega) and JM109 competent cells 
(Promega). Plasmids were extracted using Wizard Plus SV 
Minipreps DNA Purification System (Promega) and se- 
quenced by Sanger Sequencing method with the vector 
specific primers (supplementary table S1). 


Genomic Comparison of DNMT3A Locus Across 
Vertebrates 


Genomic sequences of several vertebrate DNMT3As 
[mouse (Mus musculus), human (Homo sapiens), tarnmar 
(M. eugenii), platypus (O. anatinus), chicken (G. gallus), 
and the bearded dragon (P. vitticeps)] including up to 
3,000 bases of upstream sequence of the TSS was down- 
loaded from NCBI (https:/www.ncbi.nim.nih.gov) and 
Wallabase (https:/wallabase.org). To locate conserved gen- 
omic region across vertebrates, genomic alignments and 
visualization were performed by AliTV (Ankenbrand et al. 
2017). To identify conserved transcription factor binding 
sites, the identified conserved region across mammals 
was subjected to multiple alignment by clustalW. Regions 
containing five or more highly conserved sequences around 
the TSS of DNMT3A2 across mammals were identified and 
compared with known transcription factor binding motifs 
by motif searches with JASPR CORE database (2020) using 
TomTom function in the MEME Suite (https:/meme-suite. 
org/meme/index.html; Bailey et al. 2015). Motifs with an 
E-value <0.05 and q-value <0.05 were considered as po- 
tential transcription factor binding motifs. 


Transcriptome Analysis 


Publicly available raw RNA-seq data sets of koala (GEO acces- 
sion: GSE128122; Yu et al. 2019), platypus (GEO accession: 
GSE97367; Marin et al. 2017), chicken (GEO accession: 
GSE97367; Marin et al. 2017), and the bearded dragon 
(BioProject accession: PRJEB5206; Georges et al. 2015) 
were downloaded from NCBI SRA (https:/Awww.ncbi.nim. 
nih.gov/sra). All RNA-seq reads were trimmed using 
TrimGalore! (v0.6.5; https:/github.com/FelixKrueger/ 
TrimGalore) with default settings to eliminate adaptor se- 
quences, poor quality reads, and very short (<20 bp) reads. 
The trimmed reads were aligned to each genome (koala: 
phaCin_unsw_v4.1; platypus: mOrnAnat1.p.v1; chicken: 
GRCg6a; the bearded dragon: pvi1.1) using HISAT2 (Kim 
et al. 2019). Biological replicates were merged to detect spli- 
cing variants by increasing sequence depth. The reads were 
used as inputs of ggsashimi program (Garrido-Martin et al. 
2018) to visualize splicing junctions. 


Reverse transcription Quantitative RT-PCR 


Reverse transcription quantitative RT-PCR (RT-qPCR) pri- 
mers (supplementary table S1) were designed from the 
known tammar DNMT3A and DNMT3A2 sequences and 
amplified 113 and 76 bp products, respectively. RT-qPCR 
was performed on a Quantstudio 5 (Thermo Fisher 
Scientific) on gonadal cDNA of both sexes from day 10 pp 
to adult using the SYBR Green PCR kit (QIAGEN, Hilden, 
Germany). TATA-box binding protein (TBP) and hydroxy- 
methylbilane synthase (HBS) were selected as the house- 
keeping genes. Expression levels of the DNMT3As were 
normalized to the geometric mean of the expression levels 
of housekeeping genes. Differences in the gene expression 
levels between stages were analyzed by Tukey—Kramer’s 
honestly significant difference test in R. 


Immunohistochemistry 


Extensive testing of different antibody concentrations and 
pretreatments (enzyme digestion and heat treatment in dif- 
ferent buffers) using various tammar tissues and organs 
was performed to optimize each antibody. PFA-fixed, 
paraffin-embedded wallaby gonads were sectioned at 
6 um and sections were mounted on Superfrost slides 
(Thermo Fisher Scientific, Victoria, Australia). Slides were 
deparaffinized, rehydrated, permeabilized in 0.1% (v/v) 
Triton X-100 in tris-buffered saline (TBS; 50 mM Tris HCI, 
150 mM NaCl, pH 7.4) for 15 min and then blocked in 
6% hydrogen peroxide in water for 10 min. Slides requiring 
heat treatment were heated in a water bath at 99 °C in the 
appropriate buffer (supplementary table S2) for 30 min and 
then left at room temperature for another 30 min to cool. 
Nonspecific staining was blocked using 10% serum from 
the appropriate species for 30 min at room temperature. 
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Primary antibodies were diluted as described in 
supplementary table S2 in TBS/0.1% bovine-serum albumin/ 
5% appropriate serum and then applied to the sections for 
16 h at 4 °C. Negative control sections were incubated with 
an equivalent concentration of the appropriate immuno- 
globulin. Signal was amplified using biotinylated goat anti- 
rabbit or rabbit anti-mouse secondary antibodies and then 
an ABC/HRP kit (DAKO, New South Wales, Australia) and vi- 
sualized with DAB (Sigma). Sections were lightly counter- 
Stained in hematoxylin and then coverslipped in DPX. At 
least three (usually four to five) gonads were stained for 
each stage for each antibody. Sections were viewed using 
an Olympus BX51 System microscope and photographed 
using an Olympus DP70 digital microscope camera. 


Western Blot 


Protein was extracted from day 64 pp tammar PY testis (n = 
4 pooled), day 80 pp ovary (n=5 pooled), and day 23 of 
gestation fetus (n9=2 pooled) by homogenization in RIPA 
buffer (50 mM Tris base, 150mM NaCl, 1% NP-40, 
0.25% Na-deoxycholate). Fifty micrograms of each sample 
were electrophoresed in a 7.5% SDS-PAGE gel and then 
transferred to a PVDF membrane (Amersham Hybond-P). 
The membrane was blocked in 5% skim milk dissolved in 
TBST (TBS, 50 mM Tris HCI, 150 mM NaCl, pH 7.5 with 
0.1% Tween 20) at 4 °C overnight and then incubated at 
room temperature with the anti-DNMT3A_ antibody 
(ab13888, 2 pg/ml) for 2 h, followed by a 1h incubation 
with a HRP-conjugated goat anti-mouse secondary anti- 
body diluted in TBS with 5% skim milk and 0.1% Triton 
X-100. Chemiluminescence using the Amersham ECL 
Western Blotting Detection Reagents kit (GE Healthcare 
Life Sciences) was used to detect bound antibody. 


Supplementary Material 


Supplementary figure and supplementary table are avail- 
able at Genome Biology and Evolution online (http:/ 
www.gbe.oxfordjournals.org/). 


Author Contributions 


T.I., D.H., J.C.F., O.W.G., S.S., and M.B.R. discussed and de- 
signed the research; T.I. performed most of the research; 
D.H. performed IHC experiment; J.C.F. and M.B.R. collected 
platypus tissues; J.C.F. extracted RNA from platypus tissues; 
O.W.G., S.S., and M.B.R. provided supervision; T.., J.C.F., 
O.W.G., S.S., and M.B.R. discussed the data and T.I., 
D.H., J.C.F., O.W.G., S.S., and M.B.R. wrote the paper. 


Acknowledgments 


The authors thank all members of the Wallaby research 
group for helpful discussions and assistance with sampling. 


They thank Professors Peter Temple-Smith, Geoff Shaw, 
and Dr Kath Handasyde for collecting platypus samples. 
This study was supported by a Melbourne Research 
Scholarship to Tl and grants from the Australian Research 
Council to MBR and Prof. Geoff Shaw. 


Data Availability 


Raw RNA-seq data sets of koala (GEO accession: 
GSE128122), platypus (GEO accession: GSE97367), chick- 
en (GEO accession: GSE97367), and the bearded dragon 
(BioProject accession: PRJEB5206) were publicly available 
on NCBI SRA (https:/www.ncbi.nim.nih.gov/sra). 


Literature Cited 


Ankenbrand MJ, Hohlfeld S, Hackl T, Forster F. 2017. AliTV—inter- 
active visualization of whole genome comparisons. Peer J 
Comput Sci. 3:e116. 

Auclair G, Weber M. 2012. Mechanisms of DNA methylation and de- 
methylation in mammals. Biochimie. 94:2202-2211. 

Bailey TL, Johnson J, Grant CE, Noble WS. 2015. The MEME suite. 
Nucleic Acids Res. 43:W39-W49. 

Bartolomei MS, Ferguson-Smith AC. 2011. Mammalian genomic im- 
printing. Cold Spring Harb Perspect Biol. 3:a002592. 

Bourc’his D, Xu G-L, Lin C-S, Bollman B, Bestor TH. 2001. Dnmt3L and 
the establishment of maternal genomic imprints. Science (1979). 
294:2536-2539. 

Chen T, Ueda Y, Xie S, Li E. 2002. A novel Dnmt3a isoform produced 
from an alternative promoter localizes to euchromatin and its ex- 
pression correlates with active de novo methylation. J Biol Chem. 
277:38746-38754. 

Cuneo R, et al. 2013. High-precision U-Pb geochronology and a new 
chronostratigraphy for the Canadon Asfalto Basin, Chubut, central 
Patagonia: Implications for terrestrial faunal and floral evolution in 
Jurassic. Gondwana Res. 24:1267-1275. 

Das R, et al. 2012. Convergent and divergent evolution of genomic im- 
printing in the marsupial Monodelphis domestica. BMC Genomics. 
13:394. 

Delgado I, et al. 2021. Control of mouse limb initiation and antero- 
posterior patterning by Meis transcription factors. Nat Commun. 
12:1. 

Ferguson-Smith AC. 2011. Genomic imprinting: The emergence of an 
epigenetic paradigm. Nat Rev Genet. 12:565-575. 

Frésard L, et al. 2014. Transcriptome-wide investigation of genomic 
imprinting in chicken. Nucleic Acids Res. 42:3768-3782. 

Frost JM, Moore GE. 2010. The importance of imprinting in the human 
placenta. PLoS Genet. 6:1-9. 

Garrido-Martin D, Palumbo E, Guigd R, Breschi A. 2018. ggsashimi: 
Sashimi plot revised for browser- and annotation-independent 
splicing visualization. PLoS Comput Biol. 14:1-6. 

Georges A, et al. 2015. High-coverage sequencing and annotated as- 
sembly of the genome of the Australian dragon lizard Pogona vit- 
ticeps. Gigascience 4:45. 

Glass JL, Fazzari MJ, Ferguson-Smith AC, Greally JM. 2009. CG di- 
nucleotide periodicities recognized by the Dnmt3a-Dnmt3L com- 
plex are distinctive at retroelements and imprinted domains. 
Mamm Genome. 20:633-643. 

Griffith OW, Brandley MC, Belov K, Thompson MB. 2016. Allelic ex- 
pression of mammalian imprinted genes in a matrotrophic lizard, 
Pseudemoia entrecasteauxii. Dev Genes Evol. 226:79-85. 


12 Genome Biol. Evol. 14(7) https://doi.org/10.1093/gbe/evac094 Advance Access publication 23 June 2022 


zz0z Ainr 0g uo ysanB Aq BGES199/P600eA2/Z/p1-/2/012/8q6/W09'dno‘s1Wepese//:sdyYy Woy pepeojumog 


Evolution of DNMT3A2 in Mammals 


GBE 


Hata K, Okano M, Lei H, Li E. 2002. Dnmt3L cooperates with the 
Dnmt3 family of de novo DNA methyltransferases to establish ma- 
ternal imprints in mice. Development 129:1983-1993. 

Ishihara T, Griffith OW, Tarulli GA, Renfree MB. 2021. Male germline 
development in the tammar wallaby, Macropus eugenii. 
Reproduction 161:333-341. 

Ishihara T, Hickford D, Shaw G, Pask AJ, Renfree MB. 2019. DNA 
methylation dynamics in the germline of the marsupial tammar 
wallaby, Macropus eugenii. DNA Res. 26:85-94. 

Joshi RS, et al. 2016. DNA methylation profiling of uniparental disomy 
subjects provides a map of parental epigenetic bias in the human 
genome. Am J Hum Genet. 99:555-566. 

Kaneda M, et al. 2004. Essential role for de novo DNA methyltransferase 

Dnmt3a in paternal and maternal imprinting. Nature 429:900-903. 

Keverne EB, Curley JP. 2008. Epigenetics, brain evolution and behav- 

iour. Front Neuroendocrinol. 29:398-412. 

Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. 2019. Graph-based 

genome alignment and genotyping with HISAT2 and 

HISAT-genotype. Nat Biotechnol. 37:907-915. 

a Salle S, et al. 2004. Windows for sex-specific methylation marked by 

DNA methyltransferase expression profiles in mouse germ cells. 

Dev Biol. 268:403-415. 

a Salle S, et al. 2007. Loss of spermatogonia and wide-spread DNA 
methylation defects in newborn male mice deficient in DNMT3L. 
BMC Dev Biol. 7:104. 

Lefebvre L, et al. 1998. Abnormal maternal behaviour and growth re- 
tardation associated with loss of the imprinted gene Mest. Nat. 
Genet. 20:163-169. 

Lefebvre L, Viville S, Barton SC, Ishino F, Surani MA. 1997. Genomic 
structure and parent-of-origin-specific methylation of Peg1. Hum 
Mol Genet. 6:1907-1915. 

Li T, et al. 2002. An imprinted PEG1/MEST antisense expressed pre- 
dominantly in human testis and in mature spermatozoa. J Biol 
Chem. 277:13518-13527. 

LiJ, et al. 2004. Imprinting of the human L3MBTL gene, a polycomb family 
member located in a region of chromosome 20 deleted in human 
myeloid malignancies. Proc Natl Acad Sci U S A. 101:7341-7346. 

Lucifero D, et al. 2007. Coordinate regulation of DNA methyltransfer- 
ase expression during oogenesis. BMC Dev Biol. 7:36. 

Lucifero D, Mann MRW, Bartolomei MS, Trasler JM. 2004. 
Gene-specific timing and epigenetic memory in oocyte imprinting. 
Hum Mol Genet. 13:839-849. 

Ma P, de Waal E, Weaver JR, Bartolomei MS, Schultz RM. 2015. A 

DNMT3A2-HDAC2 complex is essential for genomic imprinting 

and genome integrity in mouse oocytes. Cell Rep. 13:1552-1560. 

Marin R, et al. 2017. Convergent origination of a Drosophila-like dos- 

age compensation mechanism in a reptile lineage. Genome Res. 

27:1974-1987. 

Mattiske D, Shaw G, Shaw JM. 2002. Influence of donor age on devel- 
opment of gonadal tissue from pouch young of the tammar wal- 
laby, Macropus eugenii, after cryopreservation and xenografting 
into mice. Reproduction 123:143-153. 

Monk D. 2015. Germline-derived DNA methylation and early embryo 
epigenetic reprogramming: The selected survival of imprints. Int J 
Biochem Cell Biol. 67:128-138. 

Nimura K, et al. 2006. Dnmt3a2 targets endogenous Dnmt3 L to ES cell 
chromatin and induces regional DNA methylation. Genes Cells 11: 
1225-1237. 

Nishikawa K, et al. 2015. DNA methyltransferase 3a regulates osteoclast 
differentiation by coupling to an S-adenosylmethionine-producing 
metabolic pathway. Nat Med. 21:281-287. 

Okano M, Bell DW, Haber DA, Li E. 1999. DNA methyltransferases 
Dnmt3a and Dnmt3b are essential for de novo methylation and 
mammalian development. Cell 99:247-257. 


Ono R, et al. 2006. Deletion of Peg10, an imprinted gene acquired 
from a retrotransposon, causes early embryonic lethality. Nat 
Genet. 38:101-106. 

Pask AJ, et al. 2009. Analysis of the platypus genome suggests 
a transposon origin for mammalian imprinting. Genome Biol. 10: 
R1. 

Piedrahita JA. 2011. The role of imprinted genes in fetal 
growth abnormalities. Birth Defects Res A Clin Mol Teratol. 91: 
682-692. 

Plagge A, et al. 2004. The imprinted signaling protein XLas is 
required for postnatal adaptation to feeding. Nat Genet. 36: 
818-826. 

Rada-lglesias A, et al. 2008. Whole-genome maps of USF1 and USF2 
binding and histone H3 acetylation reveal new aspects of promoter 
structure and candidate genes for common human disorders. 
Genome Res. 18:380-392. 

Reik W, Walter J. 2001. Evolution of imprinting mechanisms: The bat- 
tle of the sexes begins in the zygote. Nat Genet. 27:255-256. 
Renfree MB, Ager El, Shaw G, Pask AJ. 2008. Genomic imprinting in 

marsupial placentation. Reproduction 136:523-531. 

Renfree MB, Hore TA, Shaw G, Marshall Graves JA, Pask AJ. 
2009. Evolution of genomic imprinting: Insights from 
marsupials and monotremes. Annu Rev Genomics Hum Genet. 
10:241-262. 

Riesewijk AM, Hu L, Schulz U, Tariverdian G, Ho P. 1997. Monoallelic 
expression of human PEG1/MEST is paralleled by parent-specific 
methylation in fetuses. Gene. 244:236-244. 

Sadakierska-Chudy A, Kostrzewa RM, Filip M. 2015. A comprehensive 
view of the epigenetic landscape part |: DNA methylation, passive 
and active DNA demethylation pathways and histone variants. 
Neurotox Res. 27:84—97. 

Saitou M, Kagiwada S, Kurimoto K. 2012. Epigenetic reprogramming 
in mouse pre-implantation development and primordial germ cells. 
Development 139:15-31. 

Sakai Y, Suetake |, Shinozaki F, Yamashina S, Tajima S. 2004. Co-expression 
of de novo DNA methyltransferases Dnmt3a2 and Dnmt3L in gono- 
cytes of mouse embryos. Gene Expr Patterns. 5:231-237. 

Sano S, et al. 2018. CRISPR-mediated gene editing to assess the roles 
of TET2 and DNMT3A in clonal hematopoiesis and cardiovascular 
disease. Circ Res. 123:335-341. 

Seisenberger S, et al. 2012. The dynamics of genome-wide DNA 
methylation reprogramming in mouse primordial germ cells. Mol 
Cell 48:849-862. 

Singh MK, et al. 2005. The T-box transcription factor Tbx15 is required 
for skeletal development. Mech Dev. 122:131-144. 

Smits G, et al. 2008. Conservation of the H19 noncoding RNA and 
H19-IGF2 imprinting mechanism in therians. Nat Genet. 40:971-976. 

Steimle JD, Moskowitz IP. 2017. TBX5: A key regulator of heart devel- 
opment. Curr Top Dev Biol. 122: 195-221. 

Stoger R, et al. 1993. Maternal-specific methylation of the imprinted 
mouse Igf2r locus identifies the expressed locus as carrying the im- 
printing signal. Cell. 73:61-71. 

Suetake |, Morimoto Y, Fuchikami T, Abe K, Tajima S. 2006. 
Stimulation effect of Dnmt3L on the DNA methylation activity of 
Dnmt3a2. J Biochem. 140:553-559. 

Suzuki S, et al. 2005. Genomic imprinting of IGF2, p57KIP2 and 
PEG1/MEST in a marsupial, the tammar wallaby. Mech Dev. 122: 
213-222. 

Suzuki S, et al. 2007. Retrotransposon silencing by DNA 
methylation can drive mammalian genomic imprinting. PLoS 
Genet. 3:e55. 

Suzuki S, Shaw G, Renfree MB. 2018. Identification of a novel antisense 
noncoding RNA, ALID, transcribed from the putative imprinting con- 
trol region of marsupial IGF2R. Epigenetics Chromatin 11:55. 


Genome Biol. Evol. 14(7) https://doi.org/10.1093/gbe/evac094 Advance Access publication 23 June 2022 13 


zz0z Alnr 0g uo ysanB Aq BGES199/P609eN9/Z/p1/20112/8q6/W09'dno‘s1Wepese//:sdyYy Woy pepeojumog 


Ishihara et al. 


GBE 


Thakur A, et al. 2016. Widespread recovery of methylation at gametic 
imprints in hypomethylated mouse stem cells following rescue 
with DNMT3A2. Epigenetics Chromatin 9:1-15. 

Ullmann SL, Shaw G, Alcorn GT, Renfree MB. 1997. Migration 
of primordial germ cells to the developing gonadal ridges in 
the tammar wallaby Macropus eugenii. J Reprod Fertil. 110: 
135-143. 

Watanabe K, Yokota K, Yoshida K, Matsumoto A, lwamoto S. 2019. 
A novel upstream transcription factor 1 target gene 
N4bp2I1 that regulates adipogenesis. Biochem Biophys Rep. 20: 
100676. 


Weidman JR, Dolinoy DC, Maloney KA, Cheng JF, Jirtle RL. 2006. 
Imprinting of opossum Igf2r in the absence of differential methy- 
lation and Air. Epigenetics. 1:49-54. 

Yokomine T, Hata K, Tsudzuki M, Sasaki H. 2006. Evolution of the ver- 
tebrate DNMT3 gene family: A possible link between existence of 
DNMT3BL and genomic imprinting. Cytogenet Genome Res. 113: 
75-80. 

YuT, et al. 2019. The piRNA response to retroviral invasion of the koala 
genome. Cell 179:632-643.e12. 


Associate editor: Dr. Rachel O'Neill 


14 Genome Biol. Evol. 14(7) https://doi.org/10.1093/gbe/evac094 Advance Access publication 23 June 2022 


zz0z Ainr oz uo ysanB Aq BGES199/P609eA9///p1/2/0112/8q6/W09'dno‘s1Wepese//:sdyYy Woy pepeojumog 


